Meson and baryon spectrum for QCD with two light dynamical quarks 



O 
(N 

(N 

-(— • 
Oh' 

(N 
> 

00 



o 
o 



X 



Gcorg P. Engel^, C. B. Lang^, Markus Limmcr^, Daniel Mohler^'^, and Andreas Schafer'^ 
(BGR [Bcrn-Graz-Regcnsburg] Collaboration) 

^ Institut fur Physik, FB Theoretische Physik, Universitat Graz, A -8010 Graz, Austria 

^TRIUMF, 4004 Wesbrook Mall Vancouver, BC V6T 2A3, Canada 

'^ Institut fur Theoretische Physik, Universitat Regensburg, D-93040 Regensburg, Germany 

(Dated: July 13, 2010) 

We present results of meson and baryon spectroscopy using the Chirally Improved Dirac operator 
on lattices of size 16'^ x 32 with two mass-degenerate light sea quarks. Three ensembles with pion 
masses of 322(5), 470(4) and 525(7) MeV and lattice spacings close to 0.15 fm are investigated. 
Results for ground and excited states for several channels are given, including spin two mesons 
and hadrons with strange valence quarks. The analysis of the states is done with the variational 
method, including two kinds of Gaussian sources and derivative sources. We obtain several ground 
states fairly precisely and find radial excitations in various channels. Excited baryon results seem to 
suffer from finite size effects, in particular at small pion masses. We discuss the possible appearance 
of scattering states, considering masses and eigenvectors. Partially quenched results in the scalar 
channel suggest the presence of a 2-particle state, however, in most channels we cannot identify 
them. Where available, we compare our results to results of quenched simulations using the same 
action. 
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I. INTRODUCTION 

The overwhelming majority of hadronic states in the 
Particle Data Group's collection are hadron excitations 
[l|. QCD as the theory of strong interactions should give 
the whole spectrum of hadrons, based on the few quark 
mass parameters and a scale. So far, the lattice regu- 
larized form of QCD provides the only known way to 
perform ab-initio calculations of these observables. Only 
recently there have been lattice calculations with dynam- 
ical quark masses close to their physical values; most cal- 
culations still rely on extrapolations from unphysically 
heavy quarks. A reliable determination of the excited 
states still remains a hard challenge. The present calcu- 
lation is another step in this enterprise. 

For decades a lattice realization of chiral symmetry has 
been an obstacle; quite early it was observed Q] that a 
solution might be Dirac operators obeying an non-linear 
relation, the so-called Ginsparg- Wilson (GW) condition, 
overlooked for many years. Meanwhile wc know several 
Dirac operators obeying that condition. One of them, 
the overlap operator [3|, |J| , has an explicit construction 
based on a domain-wall approach [5|, |6[ in the limit of 
infinite extent of a fifth dimension. Another formulation 
(perfect fermions) is formally exact [3, but, like other 
approaches, can be constructed only in an approximate 
version |8|. The lattice version of chiral symmetry un- 
derlying the GW relation has been given in Q. GW- 
Dirac operators are numerically expensive to construct 
but have nice properties like protection from additive 
mass renormalization or automatic 0{a) improvement. 

Due to the construction method, which, e.g., in the 
case of overlap fcrmion, involves computation of an op- 
erator sign function, simulations with dynamical GW- 



fermions are very expensive, typically two orders of mag- 
nitude more expensive than simulations with the simple, 
improved Wilson Dirac operator. A problem especially 
apparent for GW-fermions concerns tunneling between 
different topological sectors. Due to the Atiyah-Singer 
theorem we know that topological sectors of the gauge 
configurations are related to the net number of exact zero 
modes of the Dirac operator (counted according to their 
chirality signature). 

The so-called Chirally Improved (CI) Dirac operator 
is an app roximate solution of the Ginsparg- Wilson equa- 
tion [13, [ill ■ Its construction is based on a formal pa- 
rameterization of the Dirac operator inserted in the GW- 
equation and solved in truncated from. This fermion ac- 
tion has already been used extensively in simulations by 
the BGR-coUaboration in the quenched approximation. 
It was found that at least in quenched simulations the 
O(a^) corrections for baryon masses are small [l^ and 
that renormalization constants behave similar to an ex- 
act chirally symmetric action |13| . 

In this paper we present results of dynamical simu- 
lations with two mass-degenerate light quarks using this 
action. The parametrization, as well as details of the sim- 
ulation and some early results are given in [ij, [T5| . The 
small discretization errors allow us to use rather coarse 
lattices in order to save computational costs. 

Here we discuss the results concerning several meson 
and baryon masses derived on the gauge field config- 
urations of [13 • We present results of ground states 
as well as excited states, making use of the variational 
method [I^, [l3]- In addition to the light (dynamical 
and valence) quarks wc also consider another, heavier 
valence (strange) quark and include the strange mesons 
and baryons in our analysis. With the presently avail- 
able data we neither perform a continuum limit nor an 



infinite volume limit. This may be justified considering 
the small discretization errors and that our values for 
m,TrL should be large enough to expect small finite vol- 
ume effects due to the pion cloud. However, in addition 
to pions "wrapping around the universe" , finite volume 
effects can appear when considering large hadrons which 
may not "fit in the box". A possible interpretation of 
our results involves such squeezing effects in case of the 
excited baryons. Motivated empirically, extrapolation to 
the physical pion mass is performed naively, with a fit lin- 
ear in the pion mass squared. Preliminary results have 
been presented in |18| . 

Recent results on light and strange hadron spec- 
troscopy with focus on excited states following different 
approaches can be found in piil - l33l |. 

This paper focuses on results for hadron masses from 
dynamical Cl-simulation and is organized as follows: In 
Sect. In] we review the setup and parameters of the sim- 
ulation, followed by the discussion of the methods to ex- 
tract hadron masses in Sect. IIIII Extrapolation to phys- 
ical quark masses and sources for systematic errors are 
discussed in Sect. IIVI In Sections |V] and |Vl] we then 
present our results on hadron masses. We conclude in 
Sect. EUl 



II. SIMULATION DETAILS 
A. Fermion action and gauge action 

All details of the simulation method and our checks for 
equilibration and autocorrelation are given in [l5|. For 
completeness we summarize the essential features here. 

We use the Chirally Improved Dirac (-Dci) operator 
[lOl lll| ; this is an approximate solution to the Ginsparg- 
Wilson equation. It is obtained by insertion of the most 
general ansatz for a Dirac operator into the GW equa- 
tion and comparison of the coefficients. This leads to an 
algebraic set of coupled equations, which can be solved 
numerically. The accuracy of the approximation depends 
on the number of terms included after the truncation, in 
our case considering paths to neighbors up to a maximum 
path length of 4 lattice units. The paths and coefficients 
used are found in the appendix of [l5[. In principle, one 
could optimize the coefficients in the parametrization for 
each gauge coupling and quark mass value with respect 
to chiral symmetry. However, defining the setup this way, 
the predictive power of the simulation is weakened, and, 
furthermore, comparison of different sets of gauge ensem- 
bles is more complicated. We therefore decided to use 
the same paths and coefficients in all our dynamical runs 
and thus the bare Dirac operator is the same in all dis- 
cussed ensembles. This leads to additive mass renormal- 
ization which is corrected for by determining the PC AC 
(partially conserved axial current) mass, also called AWI 
(axial Ward identity) mass for each ensemble. 

We include one level of stout smearing [3J| as part of 
the definition of Dqi in order to improve the fermion 



action further. The parameters of stout smearing are ad- 
justed such that the value of the plaquctte is maximized 
(p ~ 0.165 in the notation of Ref. [34|'). 

It was found that the combination of the Dq\ with the 
tadpole-im pro ved Liischer-Weisz gauge action shows nice 
properties [ll[ . We use this gauge action as discussed in 



B. Algorithm 

We generate the dynamical configurations with a Hy- 
brid Monte-Carlo (HMC) algorithm Hi], with the im- 
plementation for Dq\ described in |14| . Performance im- 
provement is obtained by Hasenbusch mass precondition- 
ing [36| with two pseudofermions. Further improvement 
is achieved by the use of a chronological inverter by min- 
imal residue extrapolation [37[ and a mixed precision in- 
verter [33 • A discussion of the autocorrelation time, the 
eigenvalues of the Dirac operator and the topological sec- 
tor of the generated configurations is found in [15[ . We 
choose to analyze every fifth configuration and neglect 
any remaining weak correlations. The distribution of the 
eigenvalues of Dci indicates that we may simulate small 
pion masses on relatively coarse lattices. The algorithm 
was found to show frequent tunneling through topologi- 
cal sectors. 

The gauge configurations are determined on an SGI 
ALTIX-4700 (a machine with a peak rate of 62.3 TFlop/s 
for 9728 processors) using for each configuration (i.e., one 
unit of HMC-time) a total amount (summed over the pro- 
cessors used in parallel) of 0(60) CPUh for set A and up 
to 0(120) CPUh for set C. The quark propagators for 
the analysis are computed with a multi-mass solver and 
thus the CPU time depends only on the smallest (the dy- 
namical) quark mass. Computing the quark propagators 
for one configuration and 60 sources (five source types, 
cf.. Sect. IIIIBI and 12 Dirac-color sources each) takes 
approximately the same amount of CPU time as for gen- 
erating one gauge configuration. A recent comparison for 
different actions can be found in Ref. 1391. 



C. Simulation parameters 

We use lattices of size 16'^ x 32 at three different values 
of the gauge coupling fijj^ and bare quark mass parame- 
ter TTiQ, see Table m The lattice spacing is determined via 
the static quark potential, using a Sommer parameter of 
To,exp = 0.48 fm. The bare simulation parameters are 
chosen such that the lattice spacing is of approximately 
the same magnitude in all three ensembles. Hence their 
physical volume is of the same size as well (« 2.4 fm). 

Considering the chiral extrapolation, many calcula- 
tions make use of the so-called mass independent scheme 
(cf., the discussion in [40|)- In this scheme the lattice 
spacing is determined (for fixed bare gauge coupling) in 
the chiral limit and assigned to all ensembles with that 



set I3lw rno configs a[fni] m^[MeV] mAH'/[MeV] m^L 



A 4.70 -0.050 
B 4.65 -0.060 
C 4.58 -0.077 



100 
200 
200 



0.151(2) 
0.150(1) 
0.144(1) 



525(7) 
470(4) 
322(5) 



43.0(4) 6.4 

35.1(2) 5.7 
15.0(4) 3.7 



TABLE I: Bare parameters of the simulation: Three ensembles (A,B,C), at different gauge coupling Plw and quark mass 
parameter mo. The number of configurations, lattice spacing from the static potential assuming a Sommer parameter of 0.48 
fm, the pion mass, the (non-renormalized) AWI-mass and the dimensionless product of the pion mass with the physical lattice 
size are given. For more details see [la |. 



gauge coupling. Since so far we only have one mass value 
at each gauge coupling, we make use of a mass dependent 
scheme, differing by 0{a) corrections. 



III. METHODS IN HADRON SPECTROSCOPY 



A. Variational method 



tors is sufficiently complete. The generalized eigenvalue 
equation 



Cit)vk = Xkit,to)C{to)vk, 
Afc(i,io) oc e-(*-*°)^'' (l + 0(e-(*-*°)^^'=)) , (2) 



Excited state contributions are suppressed by 
0{c~'^^) in the hadron correlators. However, the 
so-called variational method [1^, [l3 allows to extract 
excited states in principle. This method has been used 
extensively by the BGR-collaboration [1^, l4l| - t45| and 
has gained popularity in recent years. 

A simple hadron interpolating field operator with the 
right quantum numbers will have a correlation function 
that asymptotically decays with e~^* where E denotes 
the ground state energy. However, at finite time distance 
there will be contributions from excited states embedded 
in the continuum of scattering states. On a lattice with 
finite spatial extent the corresponding energy spectrum 
is discrete. In case there are no dynamical quarks, the 
eigenstates may be interpreted as bound states of the 
valence quarks in the interpolator. In the fully dynami- 
cal situation such a simple interpretation is not possible, 
since mixing with all many particle states with the same 
quantum numbers can occur. The eigenenergy levels are 
related to the scattering phase shifts of the coupling chan- 
nels J46l - l49j and the space of lattice hadron interpolators 
has to be large enough, in order to represent the possible 
intermediate states with sufficient quality. 

Given a set of N interpolators (with given quantum 
numbers) in the variational method one constructs the 
corresponding correlation matrix 



N 

= 5^(0|O.|fc)(fc|O]|0)e- 



tEk 



(1) 



fe=l 



The idea is to offer a basis of convenient interpolators, 
wherefrom the system chooses the linear combinations 
closest to the physical eigenstates. Diagonalization of 
the correlation matrix of the interpolators thus allows to 
disentangle the physical states, if the state of interpola- 



gives the energies of the eigenstates, where AEk is the 
distance of Ek to the closest state. In the interval to < 
t < 2to it would be determined by the distance to the 
first neglected state E^+i [SO]- However, we use small 
values for to (1 or 2), since otherwise the quality of the 
diagonalization decreases. We therefore determine the 
eigenvalues in a larger window of t- values. 

The corresponding eigenvectors represent the linear 
combinations of the given interpolators which are clos- 
est to the considered physical states at each time slice. 
Hence they may be used to derive some information on 
the composition of the physical modes [5l|, |52| • 

Obviously, the number of interpolators should be large 
enough, they should be independent and have overlap 
primarily with the low modes of the theory, in order to 
reduce contamination from highly excited states. In ac- 
tual calculations, including more interpolators unfortu- 
nately increases the statistical noise in the diagonaliza- 
tion. Thus, the optimal choice is usually to use only 
those interpolators, which show good overlap with the 
low physical modes. 



B. Quark source smearing 

Hadron correlation functions are built from quark 
propagators D~^, which are computed by inverting the 
Dirac operator on a given quark source. Extended 
sources improve the signal and also allow for a larger op- 
erator basis in the variational method. We use three dif- 
ferent kinds of sources: narrow (0.27 fm), wide (0.55 fm) 
and a (P wave like) derivative source. 

The sources are computed using Jacobi smearing |53l . 
|54{ : A point-like source is smeared out by applying a 



polynomial of the hopping term, 
K 



n=0 
3 



(3) 
H{n,m) = ^{Uj{n,0)5{n + },m\ (4) 

+ Uj [n- j ,Qj 5[n- j, rhj j , 

where Sq denotes the point source. The resulting source 
shape is approximately Gaussian. The parameters k and 
K are tuned (for each ensemble of configurations) such 
as to ensure approximately the same source width in all 
ensembles. Narrow (wide) sources will be denoted by 
quark subscripts n (w) in the remainder of this paper. 

The derivative sources, Sq- , are obtained by applying 
the covariant difference operators on the wide source, S^j 

m, 

P^{x,y) = U.,{x,0)S{x + i,y)-U^{x~i,0)'<5{x~ly) , 

(5) 



'5'a 



^i^u 



where i is one of the spatial directions. The derivative 
sources were found to be crucial for some states |55| . as 
will be confirmed in this paper as well. In the following, 
derivative sources are indicated by the subscript di of the 
quark field. 



C. Constructing the interpolators 

As already mentioned, we construct several interpola- 
tors in each channel in order to be able to extract excited 
states using the variational method. All sources are lo- 
cated in a single time slice and built on configurations 
which have been HYP-smeared in the spatial directions 
three times J56| . The main motivation for link smearing 
is suppression of UV-fluctuations which manifest them- 
selves, e.g., in the distribution of the plaquette. Thus, 
the parameters of the spatial HYP-smearing have been 
optimized by a trade-off between a maximum average 
plaquette and a maximum of the minimum plaquette, 
partly following the arguments of [53, [Sg . We obtained 
the parameters ai = 0.8 and a2 = 0.4, where ai is the 
parameter in the last step of the smearing algorithm, 
where the center link is smeared. 

The center positions of the quark sources are shifted for 
subsequent configurations in order to decrease statistical 
correlation of the data. The interpolators at the sink are 
projected to zero momentum, thus for sufficiently many 
configurations the sum projects to propagators of zero 
momentum hadrons due to translation invariance. Tables 
of the interpolators are found in Appendix [X] 

There exists another approach for interpolator con- 
struction developed recently, called "distillation" jSOj . 
which we do not follow here. Recent results on hadron 
spectroscopy following this approach are found in |30l . l32| . 



1. Meson interpolators 

We consider isovector-mesons, thus there are no dis- 
connected diagrams. Using spatially isotropic sources (n 
and w) , the quantum numbers of an interpolator are de- 
termined by the combinations of the spinor components 
(Dirac content). This restricts the meson to just a few 
(non-exotic) channels of spin < 1. A way to enlarge 
the basis of interpolators and access higher spin states is 
given by considering the direct group product of spinor 
and spatial structure. The decomposition to the irre- 
ducible representations then leads to interpolators with 
definite quantum numbers |60l - [64| . We realize a non- 
trivial spatial structure by using the derivative sources, 
which transform according to the lattice spin irreducible 
representation Ti. 

Depending on the quark content and the implemen- 
tation of the derivative sources, symmetrization of the 
interpolators is needed in order to have definite C-parity. 
Hence, light meson interpolators are symmetrized prop- 
erly, while symmetrization is omitted in the strange me- 
son sector (see Appendix |X| . Our strange meson cor- 
relator calculation omits cross correlation matrix ele- 
ments corresponding to interpolators with different C- 
parity quantum numbers in the limit of degenerate quark 
masses. Therefore, when analyzing strange mesons, we 
have to restrict ourselves to subsets of interpolators shar- 
ing the same J^*^ quantum numbers in the limit of de- 
generate quark masses. 



2. Baryon interpolators 

For the construction of baryon interpolators we use 
only Gaussian smeared quark sources (n, w). In case of 
the nucleon, E and S we use three different Dirac struc- 
tures, in case of the A and Q, only one. Since a baryon is 
built from three valence quarks, there are 2'^ = 8 possi- 
ble smearing combinations. If there is isospin symmetry, 
some of the resulting 8 interpolators are very similar to 
others, which we thus prune from the considered set of 
interpolators. We end up with 18 interpolators in the 
nucleon channel, 6 in the A and f2 channels and 24 in 
the E and S channels (see Appendix \K\ . We project to 
definite parity in each channel. 



3. Energy levels 

In full QCD calculations the single hadron states cou- 
ple to channels with two or more hadrons, like the even 
number of pious in the p sector. Although the original 
hadron is projected to its rest frame, the scattering states 
have internal relative momenta. For finite spatial exten- 
sion the admissible values of the momenta depend on the 
spatial size and the (Euclidean) discrete energy levels are 
related to the phase shift of the scattering states. In the 
elastic region this relationship has been discussed in (4f 



Neglecting further interactions of the hadronic bound 
states, the energy level E{A, B) for two free hadrons 
reads 



E{A{p),B{-p)) 



(6) 



^A + W + \/< + W (1 + 0{ap)) 



The hadrons A and B have back-to-back momenta since 
the whole state is projected to zero momentum. In the 
infinite volume limit, there is a continuum of scattering 
states. In a finite box, the momentum p can take only 
discrete values, determined by the boundary conditions, 
ap = 27r(nx, Uy, Uz)/ L, where L = 16 in the present work. 
In the S wave, the lowest 2-particle state level thus shows 
vanishing relative momentum, while in the P wave, the 
lowest 2-particle state level has a momentum of a|p| = 
27r/L. 

The lowest energy levels of the fictitious two-free- 
hadron-state for each ensemble are indicated in the fig- 
ures using symbols x and -I-, provided the energy levels 
are in the range of our investigation and can be estimated 
from our results. The corresponding non-correlated sta- 
tistical uncertainty, neglecting the hadronic interaction 
and finite volume effects, is of the magnitude of 5 to 60 
MeV. For clarity, these error bars are suppressed in the 
figures if the error is smaller than 30 MeV, and further- 
more, we always omit the error bar in case of continuous 
curves of many-particle states in the figures. 

In many of our hadron correlators, due to the parame- 
ters of the simulation and the resulting pion masses, two- 
particle intermediate states will have an energy larger 
than the ground state energy. As an example, the p could 
formally couple to two pions with relative momentum (to 
build the P wave) but this is forbidden for kinematical 
reason in our case, even for the lightest pion mass of en- 
semble C. There may be a slight shift of the lowest energy 
level due to avoided level crossings, though. However, the 
higher levels could well be due to two-particle states. In 
the p sector one would expect such an energy level be- 
tween the ground state and a possible p' . In most other 
cases the expected two-particle state levels are at least 
close to other observed possible resonance state levels. 
Except for partially quenched results in the light scalar 
channel, we do not observe such scattering state levels, 
or, if they are there, we cannot distinguish them from 
single particle states. 

A possibility to shed some light on the nature of the 
state is to monitor the eigenvectors Vk of Eq. [5] of the 
state when varying parameters of the simulation. Ide- 
ally, one compares the eigenvectors for several dynamical 
simulations, but also partially quenched data can yield 
some information. Since effects from partial quenching 
can shift the energy level, corresponding results may also 
allow to extract further information about the state. 

For example, in the light and strange scalar chan- 
nels we find that our partially quenched results are well 
described by partially quenched formulae of 2-particle 



states (65|. However, at the dynamical point our data 
do not allow for a unique interpretation. Both, the reso- 
nance and the scattering state may be present and con- 
tribute to the measured energy level. 

In the case of the negative parity baryon channels, pos- 
itive parity baryons and pseudoscalars might form scat- 
tering states whose energy levels are consistent with our 
results. While the extracted masses slightly favor the 
scattering states, the eigenvectors do not allow for an in- 
terpretation in terms of a level crossing and thus do not 
confirm the picture of a scattering state at small pion 
masses, either. 

An explanation for missing scattering states would be 
weak coupling to the interpolators considered. In case of 
the S wave, there is a noteworthy amplitude already at 
small momenta, while in case of the P wave higher mo- 
menta are needed. This may explain why we see possibly 
a 2-particle signal in the S wave of oq ( J^*^ = 0++), but 
not, e.g., an additional energy level between the ground 
state and p' in the P wave of p (1 ). Consideration of 
explicit two-hadron interpolators may help. 

It is known that channels with two or more particles are 
suppressed by factors 0(1 /L^) [6^,[63]. This suppression 
comes on top of the generic suppression of the excited 
states. Even in kinematical situations where already the 
ground state couples to scattering states, it turned out, 
that one has to include both, one- and two-particle states, 
in the set of interpolators in order to see clear signals 
j68| . For such attempts, see J69l - l73l |. Unfortunately, the 
computation of cross-correlations of one- and two-particle 
states intrinsically includes disconnected diagrams, which 
are technically demanding and thus not considered here. 



IV. COMPARING TO EXPERIMENT 
A. Chiral Extrapolation 

We simulate at pion masses larger than the physical 
one. Therefore we have to rely on some extrapolation to- 
wards small quark masses, if we want to make predictions 
at the physical point. The analytic form of the chiral ex- 
trapolation depends on the path taken in the parameter 
space of, e.g., the lattice spacing, the quark mass, the 
gauge coupling and the volume. In chiral perturbation 
theory (ChPT) [tI, UM ^^e quark mass is the only vary- 
ing parameter. The lattice spacing and the volume are 
assumed to remain constant. In the mass- independent 
scheme (cf., the discussions in [40, UM) ^he scale (the 
physical value of the lattice constant) is set by extrapo- 
lating the lattice spacing towards the chiral limit along 
some path in parameter space (e.g., constant bare gauge 
coupling). The extrapolated value is then used for all 
ensembles along that path. 

Since we have only one ensemble at each value of the 
coupling, we cannot use this mass-independent scheme. 
Hence, our path of extrapolation in parameter space im- 
plies a mass-dependent scheme and formulae of ChPT 



would have to be adjusted. Nevertheless, we assume 
our path to be close to the one in the mass-independent 
scheme, and expect that the analytic form of the chiral 
extrapolation should be similar, although with different 
expansion coefficients. Therefore, we perform chiral fits 
linear in the pion mass squared for all results and discuss 
possible other fit forms in certain cases. Note that the 
fits include only the three dynamical points, the partially 
quenched points are left out. In the figures, the solid 
black curve shows the chiral extrapolation, the dashed 
lines delimit the region of one standard deviation. The 
results of the chiral fits are summarized in Fig. B^ 



B. Systematic Effects 

We set the scale via the static potential and the Som- 
mer parameter (Sect. IIICp with vq = 0.48 fm. In the 
literature also other values have been used. We have 
only one dynamical quark mass for each value of the 
gauge coupling and thus cannot extrapolate to the chiral 
limit in order to define the scale in the mass independent 
scheme. More dynamical points would be desirable for a 
more reliable chiral extrapolation and will be included in 
future calculations. 

The strange quark is considered as a valence quark 
only. In view of results including full strange quark dy- 
namics (e.g., [231) we find, at least for the ground states, 
no noticeable difference in the mass range considered 
here. 

Discretization effects have been discussed for baryons 
in the quenched simulation, where, due to the improve- 
ment of the used action, only small 0{a^) corrections 
have been identified [l4|. In oder to confirm this for 
the dynamical simulation we would have to perform our 
study at several lattice spacings and lattice volumes. 
This is a future task. In the extrapolation to the physi- 
cal quark masses, using the mass-dependent scheme, we 
disregard discretization errors and discuss finite volume 
effects only qualitatively. 

The physical size of our spatial lattice volume is w 
2.4 fm. Finite volume effects due to "pious wrapping 
around the world" are expected to be small in ensemble 
A and B (to^L > 4), however, they could be significant 
in ensemble C (tottL sa 3.7). Squeezing of large hadrons, 
which do not "fit in the box" , may be an obstacle in the 
case of excited baryons. Indeed, we seem to observe such 
effects. Studies with a larger volume are planned. 

A subtle issue is the choice of interpolators for the var- 
ious hadron states. For sufficiently large and complete 
sets and exhaustive statistics the variational method 
would produce the cigenmodcs sufficiently well. In our 
case wc arc restricted to few interpolators and modest 
statistics. We attempt to optimize this situation by 
choosing suitable subsets of interpolators, as discussed, 
motivated by plateaus in the eigenvector components and 
exponential decay of the eigenvalues in some fit range. 
This brings in certain systematic effects, which will be 



reduced only by enhancing both, statistics and the set 
of interpolators. We discuss our choice for all hadrons 
considered subsequently. Indeed, in some cases we find 
sizeable dependence on the chosen sets. 

Another possible systematic infiuence comes from 
choosing to in the variational method and the fit range 
for the generalized eigenvalues. In principle, that impact 
can be estimated by choosing several values of tg Stiid 
varying the fit range. For the final fit one should then 
choose a window where this impact is negligible. How- 
ever, in practice the corresponding choices arc restricted 
by the given signal-to-noise ratio for coarse lattices and 
weak signals. We use io = 1 in most (and to = 2 in a few) 
cases and perform a fit in the maximum range possible, 
which is, however, often limited to 3 or 4 points for noisy 
observables. 



V. RESULTS: LIGHT HADRONS 

Following the idea of the variational method, a large 
basis of interpolators should be optimal. In practice, one 
finds that increasing the number of interpolators also im- 
plies an increased statistical noise. One thus has to find 
an optimal balance when choosing the basis for the vari- 
ational method. 

In each channel we take the subset of interpolators, 
which yields the optimum plateaus of effective masses. 
In order to find this subset, we first look at the diagonal 
elements of the correlation matrix, which represent the 
autocorrelation of the interpolators used. From that data 
one has a first hint which interpolators are candidates for 
such an optimal subset. To be more concrete, we extract 
information about the interpolators such as to which 
physical state they couple dominantly, how strongly they 
are affected by contamination of excited states and un- 
til which time separation the signal of the corresponding 
physical state is reliable. 

Usually, the so obtained set of interpolators includes a 
number of interpolators which are known to have rather 
large overlap with each-another, due to, e.g., their com- 
mon Dirac structure. This in turn means that they are 
far from being orthogonal and thus not well-suited to 
build a basis of the variational method. Thus, we pro- 
ceed by singling out subsets, applying the criterion of 
maximum orthogonality, i.e., least overlap. Doing so, we 
obtain a number of candidate subsets, all of which are 
analyzed with the variational method. The fit range of 
the plateau is chosen from the first point without con- 
taminations from higher states until the plateau breaks 
down or the noise gets too large. This qualitative crite- 
rion is made more precise by choosing an optimal x^ ^-s 
quantitative criterion for the fit. 

The actual values of the energy levels are then de- 
termined by an overall exponential, correlated fit to the 
eigenvalues at all time slices in the fit range. 

This range is validated by comparison with the cor- 
responding plateau of the eigenvectors. The variational 



method simultaneously yields results for several states 
using one set of interpolators. However, we find that in 
many cases the signal can be improved by considering 
different sets for different states. The extracted masses 
of such different sets agree well within error bars, but the 
noise is reduced. Theoretically, one would expect to find 
improvement by joining both sets of interpolators, but 
in practice this means an increase of the noise, which in 
most cases is the stronger effect. We show a plot of the ef- 
fective mass of the p meson as one example for consistent 
results for two sets (see Fig. |4|). 

We present results on three ensembles of gauge con- 
figurations (cf., Table H]). In all plots, filled symbols de- 
note dynamical results and open symbols denote partially 
quenched results, where the valence quark mass is always 
larger than the sea quark mass. The symbols x and 4- 
represent energy levels of free scattering states, neglect- 
ing hadronic interaction and finite volume effects. The 
scale of the vertical axis is set by the lattice spacing, as 
discussed, unless stated otherwise. The scale of the ab- 
scissa is set using the results for the pion mass squared. 
The tables of the corresponding interpolators are found 
in the appendix. The chosen subset of interpolators is 
stated in the caption of the figures. 
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FIG. 1: Mass plot for the 0"'*' channel (pion), ground state 
and first excitation. Filled symbols denote dynamical results, 
open symbols denote partially quenched results. Parameters 
of the three ensembles A, B and C are found in Table |T] The 
list of the interpolators is found in the appendix. Results for 
the pion ground state are used to set the scale on the ab- 
scissa, here and in other figures as well. The ground state is 
measured using interpolator (1), the excitation using interpo- 
lators (3,5,9,10) in ensemble A, (1,6,9,10) in B and (3,7,8,11) 
inC. 



A. Mesons with light quarks 



2. The r 



channel: 



We simulate two mass-degenerate light quarks and 
thus use the symmetrized meson interpolators found in 
Appendix |^ 



1. The ^ channel: vr 

In Fig. [1] the mass of the pion ground state is shown as 
function of Af^ to indicate the statistical errors. In the 
variational method, the observation range of the excited 
states in this channel is limited by the backwards-running 
pion (see also the discussion in [la])- This leads to short 
plateaus for small pion masses and weakens the signal 
in our simulations. Nevertheless, we find a clear signal 
of the first excited state, compatible with experimental 
data. The signal weakens towards smaller quark masses, 
thus the fit at the dynamical point is motivated by the 
plateau at partially quenched points and the behavior of 
the eigenvalues (see Fig. ^. We find no indication of an 
intermediate three pion state signal; a Trp signal would be 
expected at even higher energy values. The excited state 
signal is expected to improve when using a lattice with a 
larger time extent, in particular in the case of ensemble 
C. We do not see a reliable and stable signal of the second 
excitation. The quenched results using the same action 
l45i showed, within the errors, similar behavior as set A. 



We find an excellent plateau for the ground state and 
an excited state signal compatible with experimental 
data (see Fig. |3]). Here, the excitation signal improved 
using io = 2, which may indicate more contamination 
from excited states in this particular channel. We de- 
cided to extract the ground state result from the set of 
interpolators (1,4), which displays a better plateau than 
the combination used for the excited state. The consis- 
tency of the ground state from (1,4) and the sets chosen 
for the excitation is shown in Fig. 01 In ensemble B we 
furthermore see a signal of the second excitation compati- 
ble with the p(1700). The physical p meson is a resonance 
which can decay into two pions with relative momentum. 
The energy of the corresponding lowest scattering state 
is determined by the mass of the pion and the minimum 
non-vanishing momentum, defined by the finite spatial 
extent of the lattice. On the lattice, for our ensemble 
parameters the energy of this scattering state would be 
above the mass of the p meson. Hence the ground state of 
the channel is dominated by the p meson, which docs not 
decay and is therefore called "stabilized" . One might ex- 
pect, however, a two-pion intermediate state between the 
p ground state energy and the signal associated to p' , but 
no such state is observed here. It is suggested that the 
coupling of the used interpolators to two-particle states 
is strongly suppressed in the P wave. 

In principle, Ti interpolators may couple to spin 3 
states with energy levels close to p' . However, from the 
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FIG. 2: Effective mass plots (Ihs) and normalized eigenvalues (rhs, logarithmic scale) for the '' channel (pion) at the 
dynamical point (bottom) and at a partially quenched point (top) of ensemble C, using interpolators (3,7,8,11). Ground state 
and the first excitation data are shown. The corresponding valence quark mass parameter and the fictitious pion mass are 
indicated in the figure. The eigenvalues suggest that the first excitation can be fitted from f = 3 to f = 7, at i = 8 we find a 
kink at partially quenched data and even a loss of the signal at the dynamical point. The effective masses are proportional to 
the derivative of the logarithm of the eigenvalues, thus showing huge error bars and very bad plateaus in case the eigenvalues 
are tumbling. The fit to the excitation at the dynamical point is motivated by the plateau at partially quenched points and 
the almost stable behavior of the eigenvalues. We remark that the final fit is an overall exponential fit of the eigenvalues for 
all points in the chosen range 3 < f < 7. 



naive continuum limit of our interpolators we expect such 
coupling to be small. 

Using only "gaussian-type" interpolators (i.e., without 
derivative sources), the energy levels of the first excita- 
tion are found to deviate slightly from the results pre- 
sented here. A possible reason could be mixing with the 
nearby higher excited states, or an early breakdown of 
the plateau which can complicate the identification of a 
correct fit range. Moving the used fit range 3-6 to 2-5 
the excitation level increases. 

Comparing with quenched results using the same ac- 
tion [45|, we find that the dynamical p ground state 
comes out significantly lighter than its quenched coun- 
terpart, which, however, is partially due to the different 
Sommer parameter value used in the quenched analysis 
(ro = 0.5 fm). Again, the first excitation of the quenched 
simulation is compatible with set A of the dynamical 
case. The dynamical points B and C indicate a steeper 
slope pointing towards the experimental results. 

The isospin singlet vector meson has mainly ss 



content and disconnected parts are suppressed ( "Zweig- 
forbidden", cf., the decay channels in the experiments). 
Thus we extract the meson mass considering partially 
quenched results of the p correlator. This is discussed in 

Sect, rvra 



3. The 0"^"*" channel: ao 

The scalar channel is difficult as usual. The mass 
plateau is short and varies depending on the chosen set of 
interpolators. The ground state mass of oq is very close 
to the lowest energy level of the (dynamica.1) two-particle 
state 7r772 (see Fig. [5]), as discussed in [78Ll79|. 

We find large effects due to partial quenching close 
to the dynamical point at small pion masses which are 
very obvious in ensemble C (see Fig. [5]). The par- 
tially quenched data do not smoothly extrapolate to 
the dy namical point. An explanation has been offered 
in [65|: The partially quenched states may couple to 
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FIG. 3: Mass plot for the 1~~ channel (p), ground state and 
first excitation. The ground state is measured using inter- 
polators (1,4), the excitation using interpolators (7,8,13,14) 
in ensemble A, (6,11,14) in B and C. The estimated energy 
level of the P wave scattering state tttt lies below the first 
excitation in all three ensembles. The scattering state is not 
observed, the reason may be too weak coupling to the used in- 
terpolators. The coupling may be especially weak in case of P 
wave scattering states. The statistical error of the 2-particle 
state (based on the errors of the particle masses involved) is 
of the magnitude of 5-10 MeV and therefore not visible in the 
figure. 




M^" [GeV ] 

FIG. 4: Zoom in the mass plot for the 1"^^ channel (p), 
ground state. The different sets of used interpolators are 
checked for consistency in the ground state. Top: Interpola- 
tors (7,8,13,14) in ensemble A, (6,11,14) in B and G. Bottom: 
Interpolators (1,4). 



pairs of pscudoscalars (composed of valence and sea 
quarks), leading to unphysical contributions that cancel 
only in the fully dynamical case. We find that our par- 
tially quenched data are well described by the partially 
quenched formulae of the scattering state, and thus inter- 
pret the partially quenched data as the 2-particle state 
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FIG. 5; Mass plot for the 0"*"^ channel. Interpolator (8) is 
used throughout. Top: Ground state of ao and chiral extrap- 
olation relying on a 1-particle interpretation are shown. Bot- 
tom: Estimate of the two-particle state nri2 for various param- 
eters shown. The blue, red and black curve (online version) 
show a prediction of the partially quenched ( "pq" ) scattering 
state Trrj2 for m^ai ^ rUsea in ensemble A, B and G, respec- 
tively. The green curve (online version) shows an estimate of 
the dynamical ("dyn") scattering state ■Kri2 (m^ai ~ rUsea)- 
Interaction of the bound states and finite volume effects are 
neglected in the energy level of the two-particle state. For 
clarity, the corresponding statistical error is omitted in the 
figure. The partially quenched data suggests an interpreta- 
tion in terms of the two-particle state 7^72, while no clear 
statement about the particle content can be made at the dy- 
namical point. 



71772. We stress that our results do not allow for a clear 
interpretation as 1- or 2-particlc state at the physical 
point. The question of the coupling of our interpolators 
to the scattering states has already been discussed in 
Sect, line 31 Note that the scattering state KK cannot 
appear in our simulations since it involves strange quark 
loops. 

In quenched simulations scattering states cannot con- 
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FIG. 6: Mass plot for the l^"*' channel (oi), ground state. 
Results using interpolators (9,10) are shown in the top plot, 
results using interpolator (8) at the bottom. The uncertainty 
stemming from the choice of interpolators is discussed in the 
text. Note that the ai(1260) is a very broad resonance. The 
estimated energy level of the 5* wave scattering state 7rp is 
below the ground state in ensemble C and comparable to the 
ground state mass in ensemble A and B. For better identifi- 
cation, we display the scattering states slightly shifted to the 
left. We may interpret that we do not observe the scattering 
state, at least in ensemble C. 



tribute to the signal when using 1-particle interpolators. 
Also, in particular in the light scalar channel, ghosts may 
appear and complicate the spectroscopy. A strategy to 
disentangle the contributions has been discussed in j4a |. 
The ground state en ergy level in quenched simulations 
with the same action [4j, |45| was extracted only at larger 
pion masses. There it was essentially compatible with our 
dynamical data of set A, extrapolating to the ao(1450) 
rather than to ao(980). The spectroscopy of the light 
scalar channel appears to benefit from sea quarks. 



^. The t^ 



channel: 



ffli 



Considering the pseudo vectors, we encounter some 
practical difficulties. Gaussian shaped interpolators do 
not yield reliable signals and derivative interpolators are 
needed. We obtain short plateaus even for the ground 
state and there appears an uncertainty associated with 
different choices of interpolators. We show results from 
two sets of interpolators to illustrate this issue (see 
Fig. ini). In principle, such a situation is possible if the 
interpolators couple to different physical states or show 
different discretization effects. In the absence of such ef- 
fects, the results from different sets should agree within 
error bars and, indeed, the error bars overlap. Note that 
the ai(1260) is a very broad resonance, hence scattering 
states may complicate the spectroscopy. Based on ex- 
periments one could expect an S wave ixp energy level 



FIG. 7: Mass plot for the 1^~ channel (6i), ground state. 
Results using interpolators (4,5) are shown in the top plot, 
results using interpolator (8) at the bottom. The uncertainty 
stemming from the choice of interpolators is discussed in the 
text. The energy level of the S wave scattering state vro; is 
estimated using the approximation m^ = rup. For better 
identification, we display the scattering states slightly shifted 
to the left. The energy level is below the ground state, in 
particular in ensemble C. 



below the ai level at least for ensemble C, which we do 
not observe. 

The energy levels of quenched simulations using the 
same action [45| appear a bit high. However, they are 
compatible with the results presented in this work if 
one considers the uncertainty associated with different 
choices of interpolators. 



5. The r channel: &i 

Here the situation is similar to the 1++ channel. The 
usage of derivative interpolators is mandatory, neverthe- 
less weak plateaus are obtained and one finds an un- 
certainty which is even larger than in the ai channel. 
Again, we show results from two combinations of inter- 
polators (see Fig. [7]). Using interpolators (4,5) one gets 
a result rather far from the experimental 6i(1235), while 
the result of interpolator (8) is compatible with experi- 
ment within a large error bar. Both results are consis- 
tent within error bars. The energy level of the S wave 
scattering state ttuj is estimated using the approximation 
m^ = nip. However, this rough estimate does not allow 
for any precise statement about the particle content of 
the measured state. 



6. The g++ channel: 



a2 



In the spin 2 channels we encounter for the first time 
the situation to have orthogonal irreducible represcnta- 
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FIG. 8: Mass plot for the 2+"*" channel (a2), ground state. Lhs; Irrep T2, using interpolators (1,2,3,4) in ensemble A and B, 
(2,3) in C. Rhs: Irrep E, interpolators (1,2,5) are used throughout. The result of ensemble C likely suffers from small statistics, 
leading to a negative slope of the chiral fit. 



tions (irrcps) on the lattice, which couple to the same spin 
state in the continuum limit. Hence we are able to com- 
pare results from these different irreps. Fig. [S] shows the 
mass of 2"*"+ in irrep T2 and E. In both representations 
using only one interpolator does not yield a reliable sig- 
nal and employing the variational method is necessary. 
The resulting mass of T2 agrees with the experimental 
02(1320) within one a. The resulting mass of E, how- 
ever, comes out too high, where the reason seems to be 
a large mass of ensemble C, leading to a negative slope 
of a linear fit. 



8. The 2 + channel: 



7r2 



Here we encounter rather large statistical errors and 
additional uncertainty associated with different sets of 
interpolators (see Fig. [TUl) Nevertheless, all results are 
consistent within 2a. Due to the large error bar, we 
find the linear chiral extrapolation compatible with both, 
the 7r2(1670) and the 7r2(1880). As in the other spin 2 
channels, the results would benefit from higher statistics. 



7. The 2 channel: p2 

Unfortunately, in the irrep T2 we can extract a mass 
only in two of the 3 ensembles (see Fig. ^. Enlarged 
statistics will be necessary in ensemble A in order to 
observe a reliable mass plateau. Also in the other en- 
sembles and in the irrep E the mass plateaus are short 
and the fit ranges are partly motivated by the clearer 
plateaus of partially quenched mass results. Neverthe- 
less, the mass obtained is consistent with the experimen- 
tal p2(1940) (see Fig. ^. An estimate of the scattering 
states appearing in this channel is shown in the figure. It 
seems possible that the measured state involves contribu- 
tions from the S wave scattering state 7ra2, however, the 
rough estimate does not allow for any clear statement. 

Note that in the continuum limit the irrep T2 couples 
also to spin 3 and irrep E to spin 4 states. Hence a signal 
of p3(1690) could be seen in T2, but not in E. However, 
the interpolators used have naive continuum limit of spin 
2, which we consider to be the reason for the missing of 
a signal of p3(1690). We stress that also in this channel 
higher statistics is necessary for a reliable extrapolation 
to the physical point. 



9. Exotics 

So-called exotic states have quantum numbers which 
cannot be constructed using a naive quark model. Most 
of the known exotic particles are found above 2 GeV, 
but also some lower ones are known, e.g. 7ri(1400) and 
7ri(1600) in the 1 ^ channel. We implemented inter- 
polators with exotic quantum numbers using derivative 
sources. Unfortunately, the obtained data are very noisy 
and we found it impossible to perform a fit of the effec- 
tive masses. Since we observe a very weak signal in the 
1 *" channel, we hope that a fit can be done when larger 
statistics becomes available. 



B. Baryons with light quarks 

Our baryon interpolators are built from Gaussian 
sources only (no derivative sources) and can be found in 
Appendix [21 All interpolators are projected to definite 
parity in each channel. 
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FIG. 9: Mass plot for the 2 channel {p2), ground state. Lhs: irrep T2, interpolator (1) is used in ensemble B and C. 
Rhs: irrep E, interpolator (1) is used throughout. The energy level of the P wave scattering state ttlu is estimated using the 



approximation mo. 



It is found to be below the ground state in all three ensembles. The higher energy level of the 5* 



wave scattering state 7ra2 is in better agreement with the observed states. However, the estimate does not allow to identify the 
particle content of the measured state. 
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FIG. 10: Mass plot for the 2 ^ channel (712), ground state. Lhs: Irrep T2, interpolator (1) is used throughout. Rhs: Irrep E, 
interpolators (1,2) are used throughout. 



1. Nucleon positive parity 

The nucleon positive parity ground state can be ex- 
tracted to good precision as usual (see Fig. [TT]). A con- 
ventional fit linear in the pion mass squared yields a 
nucleon mass a bit higher than the experimental value. 
In [80, HH it was found that a simple fit linear in the 
pion mass of the nucleon positive parity ground state 
agrees well with experiment. Indeed such an extrapola- 
tion agrees well with experiment also in our case. How- 
ever, a reliable clear distinction would only be possible 
using more data from dynamical simulations and higher 
statistics. Therefore, we stick to the expectation of the 
analytic behavior being close to the pion mass squared as 



suggested from ChPT and thus we quote the correspond- 
ing fit linear in the pion mass squared in the conclusions. 

The first excitation would be compatible with the en- 
ergy level of the P wave scattering state ttN. However, 
following the arguments of Sec. IIII C 31 we believe to see 
an almost pure 1-particle state. 

The first nucleon excitation, the so-called Roper, 
comes out several hundred MeV too high in dynamical 
simulations [1^ [3l|, [32|- There are several possible ex- 
planations for this. On one hand, whereas the negative 
parity ground state baryons are orbital excitations (ac- 
cording to the quark model), the Roper is a radial ex- 
citation. Thus its size may be substantially larger than 
that of the nucleon and affected by squeezing due to the 
spatial volume size. Another reason may be a strong in- 
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FIG. 11: Mass plot for the nucleon positive parity channel. 
Top: Resuhs for ground state and two excitations. The en- 
ergy level of the P wave scattering state nN is close to the 
first excitation in all three ensembles. For clarity, we dis- 
play the scattering states slightly shifted to the left. Bottom: 
Fit to the ground state mass of the nucleon positive parity 
channel. Interpolators (1,2,13,14) are used in ensemble A, 
(1,2,4,14,15,18) in B and C. 



flucncc from the ttN channel, which may not be properly 
represented by our set of interpolators (cf., the discussion 
in Sect. IIIIC3|) . It may be necessary to explicitly include 
such meson-nucleon interpolators, which then, however, 
poses the technical challenge of backtracking quark loops. 
Comparing to the corresponding quenched simulations 
p4| , we do not observe a significant qualitative difference 
of the results in the nucleon positive parity channel. Our 
data are also in agreement with quenched and dynamical 
results from other groups (e.g., |25l. l27l. l32| ) . 



2. Nucleon negative parity 

In general, the signal is rather noisy in negative parity 
baryon channels. In the case of the nucleon, we also find 
that the backwards-running positive parity nucleon (and 
possibly also back-to-back scattering states [82[) limits 
the fit window to the interval 3 < i < 7. However, we 
even find a signal of the first excitation, which is close 
to the ground state, both being compatible with experi- 
mental data (see Fig. [T^ . We find that all our negative 
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FIG. 12: Mass plot for the nucleon negative parity channel, 
ground state and first excitation. The energy level of the S 
wave scattering state ttA'' is close to the first excitation in 
A and close to the ground state in B and G. Naively, this 
could be interpreted as a hint for a level crossing of the 1- 
particle ground state and the scattering state. Also the low 
results in case of ensemble G suggest such an interpretation. 
However, the eigenvectors contradict this picture (see Fia:. ll3|l . 
Interpolators (2,7,9) are used in ensemble A, (1,3,7,8,9) in B 
and (1,7,8,9) in C. For clarity, we display the scattering states 
slightly shifted to the left and omit the chiral extrapolation 
in the Figure. 
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FIG. 13: Normalized eigenvectors of the nucleon negative par- 
ity channel, ground state and first excitation at the dynamical 
point. To allow for a direct comparison, the same set of in- 
terpolators (1,7,8,9) is used in all three ensembles A, B and 
G. The corresponding fits of the eigenvalues are performed in 
the time range 4 < f < 7. In all three ensembles, the ground 
state is dominated by the second Dirac structure (X2), while 
the first excitation is an almost pure xi state (see Table |X)| . 
One may conclude that no level crossing of the lowest two 
states is observed for pion masses in the range between 320 
and 520 MeV. Note that in order to extract masses, different 
sets are used in A and B, which yield more stable plateaus. 
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parity baryon ground states come out somewhat too low. 

In nature, the S wave state ttN (assuming free parti- 
cles) is below the one-particle ground state in the nucleon 
negative parity channel. This may also be true at slightly 
larger pion masses, e.g., in ensemble C. That would ex- 
plain the low results in C and thus the low chiral extrap- 
olation as well. At large pion masses, the scattering state 
becomes heavier than the 1-particle ground state and one 
expects a level crossing to take place in between. Indeed, 
our results on masses in the three ensembles are compat- 
ible with such a picture. 

However, in contrast to the scalar channels, we can 
extract useful information from the eigenvectors in the 
nucleon and sigma negative parity channels. The reason 
is that several interpolators are needed for a good sig- 
nal and that two states arc observed. Wc find that in 
all three ensembles the ground state is dominated by the 
second Dirac structure (X2), while the first excitation is 
an almost pure xi state (see Fig. [13] and Table IX)) . This 
property is seen even more clearly for partially quenched 
points, where the plateau is more stable. One may con- 
clude that no level crossing of the lowest two states is 
observed for pion masses in the range of 320 to 520 MeV. 

Another hint comes from the comparison with old 
quenched results using the same action [4J]. Using only 
3-quark nucleon interpolators, no scattering states can 
appear in quenched simulations. Hence, the eigenvec- 
tors in the quenched simulations can clearly be identified 
with 1-particle states. We stress that in the quenched 
approximation ghosts appear at low pion masses, thus a 
reliable comparison to dynamical simulation can be done 
only at large pion masses. We assume the pions of en- 
semble A to be heavy enough to allow for a comparison 
with quenched simulations. By subsequent comparison of 
dynamical simulations at decreasing pion masses, in prin- 
ciple conclusions down to physical pion masses are pos- 
sible. Unfortunately, direct comparison of the quenched 
eigenvectors with our dynamical ones is impossible, since 
the corresponding interpolators differ slightly (e.g., in the 
width of the Gaussian source). However, qualitatively, we 
find the same Dirac structure content of the two lowest 
states. From that, one may conclude that both our lowest 
two states are 1-particle states, which is also suggested 
from the missing level crossing of the eigenvectors. 

So, while the extracted mass values slightly favor the 
interpretation of a scattering state in ensemble C, the 
eigenvectors tell a different story. Since the line of argu- 
ments based on the eigenvectors seems to be more reli- 
able, we believe to sec almost pure 1-particle states and 
quote the corresponding chiral extrapolation in the sum- 
mary. 

In the quenched simulation with the Cl-action (44l | , nu- 
cleon negative parity masses have been extracted only at 
larger pion masses. The bending down at low pion masses 
observed in the present work can thus not be compared 
directly. 

Note that in order to extract masses, different sets arc 
used in A and B than the ones shown in Fig. [T31 In en- 
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FIG. 14: Mass plot for the A positive parity channel, ground 
state and first excitation. The energy level of the P wave 
scattering state nN is close to the first excitation in A and 
B, but clearly separated from any measured state in C. In- 
terpolators (1,2,3) are used in ensemble A and C, (2,4,6) in 
B. 



semble C, the fit range is motivated by tracing of plateaus 
of partially quenched results. 



3. A positive parity 

In the A positive parity channel, wc find the masses 
being too high, in particular the first excitation (see 
Fig. [H)) . although they are lower than in the quenched 
analysis [4J|. Since the statistical error is fairly small, 
systematic errors and finite volume effects seem to be re- 
sponsible. As will be discussed below, this channel - for 
large values of the valence quark mass - is used to identify 
the strange quark mass parameter. Analogous to other 
P wave scattering states, ttN seems to be missing in our 
simulation. 

Quenched results using the Cl-action [4J| have shown 
a similar systematic upwards shift of the masses. Our 
A positive parity ground state is compatible with other 
groups (e.g., [S^l)- However, there the first excitation is 
fairly close to the ground state, which we do not observe 
here. A possible reason could be the larger basis of in- 
terpolator used in that work. Hence, we cannot exclude 
possible systematic errors associated with the choice of 
interpolators in this channel. 



4- A negative parity 

There is a clear signal found of the A negative parity 
channel ground state (see Fig. [TS]). The energy level of 
the S wave scattering state ttN is below the ground state 
in all three ensembles. The results seem to be in better 
agreement with an interpretation in terms of a 1-particle 



15 



3.0 



2.5 



O 
--2.0 



1.5 



1 ' 


1 






1 






1 


1 A(-) 




/ 










1351 fl^i 7r\ - 




■■'/■■'' 






B i 


EBffi) 




- 


^ 


m^^ 


i 


* 








- 


.4*y 












• 


A 


M^% 












■ 


B 


^■-0^ X 












♦ 


C 














* 


Exp 


X 












X 


A:(JtN) 














X 


B:(JlN) 














X 


C:(JlN) 


1 


1 






1 






1 



).0 0.2 



0.4 



0.6 0.8 

M ^ [GeV^] 



1.0 



1.2 



FIG. 15: Mass plot for the A negative parity channel, ground 
state. The energy level of the S wave scattering state -nN is 
below the ground state in all three ensembles. We display the 
scattering states slightly shifted to the left. Interpolator (3) 
is used in ensemble A, (2,3) in B and C. 



state. 

Quenched results of [4J| have been extracted only at 
rather large pion masses, making a comparison to present 
results uninstructive. 



VI. RESULTS: STRANGE HADRONS 

We extract the strange quark mass parameter by iden- 
tification of a partially quenched A (i.e., valence quark 
of larger mass) with i7(1670), as described in Sect. IVI Al 
The interpolators used for strange meson spectroscopy 
are listed in Appendix [Xj Since C-parity is an exact 
symmetry for mesons only in the case of mass-degenerate 
quarks, we consider interpolators without projection to 
definite C-parity for the strange mesons. This means 
that we do not perform a symmetrization of the interpo- 
lators in the strange meson sector. Our strange meson 
correlator calculation lacks cross correlations between in- 
terpolators of different C-parity quantum numbers in the 
limit of degenerate quark masses. Therefore, when an- 
alyzing strange mesons, we have to restrict ourselves to 
subsets of interpolators sharing the same J^ quantum 
numbers in the limit of degenerate quark masses. 

In all plots shown for strange hadrons, the full sym- 
bol represents a hadron where the valence strange quark 
mass is determined from the 17(1670) and the light quark 
has the mass of the dynamical quarks. The open sym- 
bols denote partially quenched data, where only the light 
valence quark mass varies. We thus neglect the effect of 
a dynamical strange quark, motivated by the dominance 
of light quark loops over strange quark loops. 

A quantum field theory, where the sea quark masses 
do not agree with the valence quark masses, is "sick", 
as can be seen, e.g., by the appearance of ghosts which 
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FIG. 16: Extracting the strange quark mass parameter by 
identifying a partially quenched A with 51(1670). The lowest 
(magenta) horizontal line represents the n(1670). Crossing 
this horizontal line with the partially quenched A mass curves 
of ensembles A, B and C defines the strange quark mass pa- 
rameter of the three ensembles, illustrated by the three verti- 
cal lines in the corresponding colors. The resulting bare quark 
mass parameters are -0.020, -0.015 and -0.022, for ensembles 
A, B and G, respectively. The calculation for the excited f2 
around 2400 MeV is expected to be too high due to finite 
volume effects. 



violate the spin- statistic theorem. As discussed in j65l |. 
the correlators are strongly affected by partial quenching 
if the valence quarks are lighter than the sea quarks. In 
our simulations we consider only valence quarks heavier 
than or equal to the sea quarks. 

Clearly, it would be desirable to include strange sea 
quarks in the simulation. 



A. Setting the strange quark mass 

The Sl(1670) consists of three strange valence quarks 
and shows weak dependence on the light quark masses. 12 
and A share the same J^ quantum numbers, they differ 
only in their flavor content. Therefore, we use our par- 
tially quenched results in the A positive parity channel 
to identify the strange quark mass parameter for our en- 
sembles A, B and C (see Fig. [11]). We decided to choose 
strange quark mass parameters which fit the experimen- 
tal n within our error bars, allowing only for parameter 
values in discrete steps of 0.05 additive to the sea quark 
mass. In case of ensemble A the obtained mass parame- 
ter perfectly fits one of the already available quark prop- 
agators. In case of ensembles B and C, we decided to 
recompute the quark propagator at the strange quark 
mass instead of interpolating between the two adjacent 
values. 

Another possibility to extract the strange quark mass 
parameter would be to apply the same recipe to a par- 
tially quenched (f> in the 1 channel. The decay channels 
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FIG. 17: Cross-check of the obtained strange quark mass pa- 
rameter: The partiaUy quenched from the ground state of 
the 1 channel fits the experimental 0(1020) very nicely. 
The result for the excited is higher than the experimental 
value, the deviation may be due to the neglected disconnected 
diagrams or simply due to the weak signal. 



of (/)(1020) suggest that it is dominated by its {ss) flavor 
content. We use as a cross-check for the strange quark 
mass obtained via ri(1670). The result fits the experi- 
mental 0(1020) very nicely (see Fig. fTT)) . indicating that 
our approach is consistent. The excited state of 51 is as- 
sumed to suffer from finite size effects. The signal of the 
excitation of suffers from neglected disconnected dia- 
grams and poor statistics. Thus these two levels are not 
appropriate for further checks of the strange quark mass. 
However, the ground state levels of S and S positive par- 
ity may be regarded as additional affirmative cross-checks 
(see Figs. [25] and [271). 

Even in the case of the fJ, one could expect finite size 
effects, since they show up in the A positive parity chan- 
nel for both, ground and excited state. However, using 
heavier quarks, finite size effects are expected to become 
less and less important in the ground state. The cross- 
check using verifies this expectation. 



B. Mesons virith strange quarks 

1. The 0~ channel: K 

In the strange meson channel 0~ we find a very accu- 
rate determination of the K ground state, and a fairly 
reliable result of the first excitation which is compatible 
with experiment (see Fig. [T5)) . Comparing to the pion 
channel, we find that the signal improves in case of a 
heavier valence quark, as expected. 

Comparing to the corresponding quenched simulations 
[431, we do not observe a qualitative difference. 



FIG. 18: Mass plot for the 0~ channel {K), ground state 
and first excitation. Ground state measured using inter- 
polator (4), excitation measured using (4,8,11,14,17) in A, 
(4,7,11,14,17) in B and C. Note that the strange meson inter- 
polators are not symmetrized (see Appendix[X]). The P wave 
scattering state ttK* is assumed to be suppressed and there- 
fore not indicated, since at the used simulation parameters its 
energy is comparable to the one of the first excitation. 



2. The 0+ channel: Kq 

The strange scalar channel is peculiar similar to the 
light scalar channel. Furthermore, the lowest experimen- 
tal state in this channel, A'q(800) or k, is not established 
and a very broad resonance with a width of more than 
80% of its mass. It is thus not clear whether this state is 
expected to be observed in our simulation. The lowest es- 
tablished resonance listed by the Particle Data Group is 
Kq{14:30). At the parameters we use, an appearance of a 
low-lying scattering state ttK is possible. The lowest mo- 
mentum two-particle state (assuming free particles) is in- 
dicated in the figure, together with the result and the ex- 
perimental K^{800) and K^{U30) (see Fig.[Tll). The re- 
sult is compatible with both, the resonance A"q(800) and 
the scattering state. Hence we cannot definitely exclude 
either of the possibilities. Comparison to the results of 
the light scalar channel 0++ (oq) may suggest that the 
scattering state contributes at least at partially quenched 
points. Higher statistics and another volume would be 
desirable in order for a more clear distinction. If we use, 
e.g., interpolators (12,13), the eigenvectors look very sim- 
ilar for different valence light quark masses, which means 
that the state we extract remains roughly the same over 
the range of partial quenching we investigate. This can 
be interpreted as an indication that the level crossing of 
a scattering state and the resonance, if it exists, is not 
located in the mass range of our investigation. 
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FIG. 19: Mass plot for the 0"*" channel, measured with inter- 
polator 13. Since ^^0(800) is a very broad resonance, thus 
we show its experimental width only. Top: Ground state 
and chiral extrapolation according to a 1-particle interpre- 
tation (Kq) are shown. Bottom: The crosses represent an 
estimate of the dynamical S wave 2-particle state nK. The 
blue, red and black curve (online version) show a prediction 
of the corresponding partially quenched ( "pq" ) state nK for 
m-vai ^ rnsea in ensemble A, B and G, respectively. Since 
there are no strong effects from partial quenching in this chan- 
nel, the partially quenched prediction almost hits the dynam- 
ical one. Note that in the partially quenched case only the 
light valence quark varies and the strange valence quark mass 
is held fixed. Thus, considering the scattering state, the par- 
tially quenched pion mass varies, while the K mass is constant 
along the partial quenching. The results do not allow for a 
clear statement about the particle content of the state. 



3. The 1 channel: K' 

We obtain an accurate result of the ground state 
iiT* (892) and a signal of an excited state (see Fig. [^ . 
Interestingly, the excitation fits the iiT* (1680) way bet- 
ter than the ii'*(1410). Compared to the light meson 
sector, mixing of J^"*" and J^~ interpolators can ap- 



2.5 



2.0 



O 
g 1.0 



0.5 




'•%o 



I , I 

0.2 0.4 

M ' [Gev'l 



A 
B 
C 
Exp 

A:(JiK) 
B:(jiK) 
C:(JlK) 

\ 

0.6 



FIG. 20: Mass plot for the 1~ channel {K*), ground state and 
excitation. Ground state measured using interpolators (1,5), 
excitation measured using (1,5,17,18,19). The energy level of 
the P wave scattering state ttK is well separated from any 
observed state in all three ensembles. 



pear in the strange meson sector, which is expected to 
shift the energy levels. This effect is due to the miss- 
ing isospin symmetry when the involved quarks are not 
mass-degenerate and thus it grows when increasing the 
mass difference between the light and the strange quark. 
However, 1 ^ interpolators are exotic, hence it is not 
clear how much they contribute to the final state. It 
is possible that our interpolators do not show overlap 
with the ii'*(1410), because wc restrict ourselves to the 
non-symmetrized interpolators corresponding to the light 
meson channel 1 . Another possible reason would be a 
tetraquark dominance of the i^*(1410) (which is gener- 
ally not expected) or just too small statistics. However, 
quenched results using bilinear quark sources have been 
compatible with the i^*(1410) ^. 

We are planning to investigate this channel more thor- 
oughly using a larger basis in the variational method and 
a larger set of ensembles in the future. We do not observe 
the P wave scattering state 7r/\ in our simulation. 



4-. The 1^ channel: Ki 

Similar to the strange 1^ channel, we seem to miss the 
ground state i^i(1270) in the strange 1"*" channel (see 
Fig. [2T|) . Our result is in good agreement with the ex- 
perimental first excitation 7^i(1400). In this channel the 
situation is clearer, due to missing isospin symmetry, one 
expects a mixing from the l"*""*" and 1"' interpolators. In 
contrast to the latter case, these are both non-exotic. So 
far, we cannot say how strong the mixing is in our sim- 
ulations, however, it was found that the mixing is weak 
down to pion masses of 400 MeV [30|- Again, a more 
thorough investigation is hoped to shed some light on 
this issue. The results of ensemble C suggest that we do 
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FIG. 21: Mass plot for the 1"'" channel (Ki). Measured using 
interpolators (13,14,15). The energy level of the S wave scat- 
tering state nK* is close to the ground states of A and B, but 
much lower than the ground state of C. 



not observe the S wave scattering state ttK* 



5. The 2 channel: K2 

As already discussed, there are two different orthogo- 
nal lattice irreducible representations for spin 2 channels. 
Again, we analyze the irrcps independently and discuss 
the results. Similar to the corresponding light meson 
channels, we encounter large error bars in the strange 
2~ channel (see Fig. E^ . The results are compatible 
with both the experimental ii'2(1770) and the K2{IS2Q) 
in both representations, albeit with an error of slightly 
more than one a in case of T2. The X2(1580) is not 
completely confirmed experimentally so far. The lowest 
possible scattering state [ttK* in P wave) is not observed 
in our simulation. Like in the spin 1 channels, we hope 
that enlarging the basis will improve the signal. 



6. The 2+ channel: Ki 

In the 2+ channel the signal is somewhat better than 
in the 2~ channel. Here we find the result of irrep T2 
closer to the experimental _fi'|(1430) than the one of ir- 
rep E (see Figs. [23l and [24| . However, the negative slope 
of the chiral extrapolation in irrep E is clearly an arti- 
fact of too small statistics (remember that only the three 
dynamical points enter the linear fit). A more complete 
basis and additional ensembles are expected to improve 
the calculation in this channel. 



7. The 1 channel: (f> 

The decay channels of suggest that it is an almost 
pure ss-state. Therefore, in the partially quenched anal- 
ysis we can read off our results for from the partially 
quenched data in the 1 [p) channel without any chi- 
ral extrapolation (see Fig. [T7)) . As already discussed, 
the ground state reproduces (j){\Q2Q) nicely, which can 
be seen as an affirmative cross-check for our method to 
set the strange quark mass parameter. Reading off the 
excited state we find some deviation from the experimen- 
tal value, explicitly shown in Fig. [29] in the summary. 
The origin of this discrepancy may be due to neglected 
disconnected diagrams or just lie in the weakness of the 
corresponding effective mass plateau. Difficulties with 
excitations are found also in related channels. E.g., in 
the 1^ channel we miss the first excitation of K* and in 
the light 1 channel the excitation of p is quite noisy. 



C. Baryons with strange quarks 

1. E positive parity 

The interpolators of S and S have the same Dirac 
structure as the nuclcon interpolator, they just differ by 
the flavor content. Hence we use similar sets of interpola- 
tors in the variational method. We obtain a ground state 
and two excited states in the E positive parity channel 
(see Fig. [5n|) . The ground state result is in satisfactory 
agreement with experimental data. This is another af- 
firmative cross-check of the setting of the strange quark 
mass parameter. The excitations are too high, which 
may be due to finite volume effects, as we have discussed 
already. Analogously to previous channels, we do not 
believe to see contributions from the P wave scattering 
state KN. 



2. E negative parity 

We find a ground state and two excitations in the E 
negative parity channel (see Fig. [521). Similar to the nu- 
cleon negative parity channel, the first excitation is very 
close to the ground state, which is in qualitative agree- 
ment with the experimental E(1620) and S(1750). The 
ground state mass agrees with experiment within error 
bars, though being a bit low, as in all other our negative 
parity baryon channels. The second excitation could be 
compatible with the observed I](2000), after correction 
of finite volume effects. I](1620), I](1750) and I](2000) 
are classified with 2, 3 and 1 stars, respectively, by the 
Particle Data Group [1|. We can confirm the two lower 
ones qualitatively and the existence of an excitation in 
the vicinity of 2000 MeV. 

In nature, the S wave scattering state KN , which is 
lighter than the 1-particle ground state, appears in the 
E negative parity channel. We analyze the eigenvec- 
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FIG. 22: Mass plot for the 2~ channel (^^2), ground state. Lhs: Irrep T2, measured using interpolators (1,3). Rhs: Irrep E, 
measured using interpolators (1,3). The experimental A'2(1580) is not completely confirmed so far. The energy levels of the 
scattering states 7rA'2 [S wave, dominant in the experiment) and nK* [P wave) are shown. 
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FIG. 23: Mass plot for the 2+ channel (JvJ) in irrep T2, 
ground state. Measured using interpolators (1,7) in A, (3,5) 
in B, (2,4,5,6) in C. 



FIG. 24: Mass plot for the 2+ channel (iv|) in irrep E, ground 
state. Measured using interpolators (2,3,9) in A and C and 
(3,4,9) in B. 



tors analogously to the nucleon negative parity channel. 
Again, wc find that in all three ensembles the ground 
state is dominated by the second Dirac structure (X2), 
while the first excitation is an almost pure xi state (see 
Table IX)) . The partially quenched points show the same 
behavior. One may conclude that no level crossing of 
the lowest two states is observed for pion masses in the 
range of 320 to 520 McV. As in the nucleon negative par- 
ity channel, wc find that the low results for masses at 
small pion masses may be explained by the presence of a 
scattering state, but the eigenvectors do not confirm this 
picture. Since the argumentation based on the eigenvec- 
tors seems to be more reliable, we believe to see an almost 
pure 1-particle state and quote the corresponding chiral 
extrapolation in the summary. 



3. H positive parity 

In the S positive parity channel we obtain a ground 
state with rather small error which is in very good agree- 
ment with the experimental S ground state (see Fig. l27p . 
We also get a prediction for a first excited state with com- 
paratively small error bar in the range of 2200 to 2400 
MeV, which after correction of finite volume effects could 
be compatible with the S(2120) (1 stars) or the S(2250) 
(2 stars), where both states are listed by the particle 
data group stating neither spin nor parity. Analogous 
to previous cases, we do not believe to see measurable 
contributions from the P wave scattering state KH. 
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FIG. 25: Mass plot for the E positive parity channel, ground 
state and two excitations. Ground state measured using inter- 
polators (1,6,17,20), first excitation measured using (1,6,20), 
second excitation measured using (3,4,8,9,13). The energy 
level of the P wave scattering state NK is very close to the 
first excitation in all three ensembles. For better identifica- 
tion, we display the scattering states slightly shifted to the 
left. 



4-. 2 negative parity 



FIG. 26: Mass plot for the E negative parity channel. Mea- 
sured using interpolators (1,9,10,12). The energy level of the 
S wave scattering state KN is very close to the ground state 
in all three ensembles (also to the first excitation in case of 
ensemble A) . For better identification, we display the scatter- 
ing states slightly shifted to the left. The results would be 
compatible with an interpretation in terms of a level crossing 
of the 1- and 2-particle states. However, analogously to the 
nucleon negative parity channel, the eigenvectors contradict 
this picture. For clarity, the chiral extrapolation is omitted 
in the figure. 



In the negative parity channel there is not even one 
established state listed by the Particle Data Group. We 
observe a ground state and a close first excitation in the 
range of 1550 to 1800 MeV, both with relatively small 
error (similar to the quenched results 14J]). These states 
may match the S(1620) (1 stars) and the S(1690) (3 
stars). There is a signal of a second excitation around 
2400 MeV as well, albeit with rather large statistical un- 
certainty and suffering from finite volume effects. We 
found it crucial to include all three baryon Dirac struc- 
tures in the variational method in order to obtain both 
excitations. 

Compared to quenched results of [4J] , we extracted an 
additional excited state in the present work. 
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5. n positive parity 

As already discussed, we use our result for fi(1672) 
to set the strange quark mass parameter (see Fig. [TO]) . 
Similar to the case oi (p, we read off results for il from 
partially quenched data for A, without any chiral extrap- 
olation. In addition to 51(1672), we obtain a first excita- 
tion in the positive parity channel (explicitcly shown in 
Fig.[29]in the summary). Since there is no extrapolation 
involved, its statistical error is rather small, however, the 
true mass value may be somewhat smaller due to finite 
volume effects. The Particle Data Group lists £7(2250) (3 
stars), fl(2380) (2 stars) and rj(2470) (2 stars), stating 
neither spin nor parity. Taking into account finite vol- 



FIG. 27: Mass plot for the H positive parity channel. Mea- 
sured using interpolators (3,4,9,13,16). The energy level of 
the P wave scattering state KT, is very close to the first ex- 
citation in all three ensembles. 



umc effects, the result for the excitation may turn out to 
be compatible with 17(2250). 



6. Q negative parity 

The Particle Data Group does not list any established 
state in the fl negative parity channel. We observe a 
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0.4 
M ^ [GeV^] 



FIG. 28: Mass plot for the H negative parity channel. 
Measured using interpolators (4,9,10,14,23) in A and B, 
(3,7,9,10,14) in C. For clarity, the chiral extrapolation is omit- 
ted in the figure. 



ground state in the range of 2050 to 2100 MeV (compare 
Fig. [151 explicitely shown in Fig. [29] in the summary). 
Such a state has also been observed in the quenched study 



44| . It does not fit to any state listed by the Particle 



Data Group. 



VII. CONCLUSION 



We have presented results of hadron spectroscopy us- 
ing the Chirally Improved Dirac operator on lattices of 
size 16^ X 32 with two mass-degenerate light sea quarks. 
Three ensembles with pion masses of 322(5), 470(4) and 
525(7) MeV and lattice spacings all close to 0.15 fm have 
been investigated. This allows for a naive chiral extrap- 
olation in the mass-dependent scheme, but neither for 
a continuum nor a thermodynamic limit. We have dis- 
cussed possible systematic effects in Sect. [TV] Systematic 
uncertainty due to discretization effects is not explored 
explicitely and may be non-negligible for some obscrv- 
ables. Finite volume effects are discussed only qualita- 
tively, where we find some indications in case of excited 
baryons. 

We have shown results for ground states and excited 
states for several meson and baryon channels, includ- 
ing spin 2 mesons constructed by the use of derivative 
sources. The spectrum of strange hadrons was accessed 
by using partially quenched strange quarks. Possible ef- 
fects from partial quenching have been discussed briefiy 
and they seem to play no important role for most of our fi- 
nal results. However, including strange sea quarks in the 
simulation would be desirable in order to reduce possible 
sources of systematic errors. The value of the strange 
quark mass parameter was set by identification of the 
partially quenched positive parity A with the 57(1672). 

The results are summarized in Fig. 1291 Several of the 
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FIG. 29: Collection of the results for light meson masses (top), 
strange meson masses (middle) and baryon masses (bottom). 
Indicated errors are of purely statistical nature. Results are 
obtained by chiral extrapolation of dynamical data linear in 
the pion mass squared (except for (j) and Q hadrons). Ex- 
perimental values listed by the Particle Data Group [J] are 
denoted by horizontal lines, the ones needing confirmation 
by dashed lines. Circles indicate ground states, squares first 
excitations and diamonds second excitations. Results shown 
aside each other are obtained using different sets of interpola- 
tors aiming for the same state. In the spin 2 channels, results 
from T2 interpolators are shown on the left, results from E 
on the right. In the p2 channel we have results only from 
E interpolators. Strange quarks are implemented by partial 
quenching, the strange quark mass parameter is set by iden- 
tification of a partially quenched A with f2(1670). Excited 
baryons which come out too high seem to systematically suf- 
fer from finite volume effects. Higher statistics will be needed 
to establish if the results for the ground states in the nega- 
tive parity baryon channels are compatible with experiment 
or whether the results are systematically too low. Results 
for scalar mesons are omitted in this figure, since the particle 
content of the corresponding observed states is too uncertain. 
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experimentally known ground states are reproduced with 
fairly high precision {p, K^ , K* , N and S positive par- 
ity). In addition, various radial excitations are found 
(tt, p, K* , A positive parity, N, E and S both pari- 
ties) . We discussed the possible appearance of scattering 
states in various channels. The coupling of our inter- 
polators to many-particle states seems to be weak and 
such states are barely, if at all, identifiable. Only in the 
light scalar channel the partially quenched data suggest 
a large contribution from an S channel 2-particle state of 
pseudoscalars. However, at the dynamical point no clear 
statement is possible. In the negative parity nucleon and 
S channels, the eigenvectors do not confirm the picture 
of the S wave 2-particle states, either, although such an 
admixture cannot be completely excluded. A clear inter- 
pretation of the particle content of the observed states is 
still missing. 

Where possible, we discussed the influence of sea 
quarks by comparison to quenched results obtained from 
the same action. In particular the spectroscopy in the 
light scalar channel seems to benefit from dynamical 
quarks, nevertheless it remains a difficult channel. Also 
the results in the light vector channel are in slightly bet- 
ter agreement with experiment when sea quarks are in- 
cluded. However, in most channels we did not observe 
a significant difference between quenched and dynamical 
simulations. We stress that in most cases comparison is 
difficult since some other details of the simulations differ 
as well. Also we have to emphasize that in many observ- 
ables the overall effect of the sea quarks may approxi- 
mately cancel (e.g., by including their effect on the scale 
setting). In other words, different dynamics may partly 
show similar phenomenology. A similar mechanism pos- 
sibly works in some of the strange hadrons, where our re- 
sults are in good agreement with experiment, despite the 
partial quenching approximation of the strange quark. 

Several channels are expected to benefit from enlarged 
statistics, especially radial and spin excitations, the pseu- 
dovector mesons and negative parity baryon channels. 
In the nucleon positive parity channel we observe excita- 
tions which arc definitely higher than the expected Roper 
resonance. This problem may be due to finite size effects 
or the lack of explicit 2-particle interpolators in the anal- 
ysis. Our data confirm the existence of some states of un- 
clear status and predict some states which are not listed 
by the Particle Data Group so far, albeit with large sta- 
tistical errors and also non-negligible systematic errors in 
some cases. 
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Appendix A: Tables of interpolators 
1. Meson interpolators 

We list the used interpolators for each meson channel 
in TablcsHIl- IIXI In case of the light mesons, symmetriza- 
tion is performed for definite C-parity. Our strange me- 
son correlation calculation lacks cross correlation ma- 
trix elements according to interpolators with different C- 
parity quantum numbers in the limit of degenerate quark 
masses. Therefore, when analyzing strange mesons, we 
have to restrict ourselves to subsets of interpolators shar- 
ing the same J^ quantum numbers in the limit of de- 
generate quark masses. 

The letters a and b denote light or strange quarks; 
n [w) denotes a Gaussian shaped narrow (wide) source 
and di denotes a derivative source in spatial direction i 
(see Sec, imp . 74 is the gamma matrix in time direction, 
tijk is the Levi-Civita symbol, Qijk are Clebsch-Gordon 



coefficients, where all elements are zero except Qm 

72' ^122 = -^ 
2 

x/6- 



^211 — ^"Te' '3222 



-^ and (5233 = 



Interpolator 
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I0++ 
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2o++ 
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^W ^\D 
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60++ 
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TABLE III: Scalar interpolators from irrep Ai. The un- 
symmetrized interpolator, the corresponding number in the 
strange meson channel 0^, the number and the symmetriza- 
tion in terms of the 0~ interpolators are given for the light 
meson channel 0^^. 
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TABLE II: Pseudoscalar interpolators from irrep Ai. The un- 
symmetrized interpolator, the corresponding number in the 
strange meson channel 0~, the number and the symmetriza- 
tion in terms of the 0~ interpolators are given for the light 
meson channel 0""^. 
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TABLE IV: Vector interpolators from irrep T\. The un- 
symmetrized interpolator, the corresponding number in the 
strange meson channel 1~, the number and the symmetriza- 
tion in terms of the 1~ interpolators are given for the light 
meson channel 1 . 
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Interpolator 
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TABLE V: Pseudovector interpolators from irrep Ti. The unsymmetrized interpolator, the corresponding number in the strange 
meson channel l""", the number and the symmetrization in terms of the 1^ interpolators are given for the light meson channels 
1++ and 1+". 
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TABLE VI: Pseudotensor interpolators from irrep T2. The unsymmetrized interpolator, the corresponding number in the 
strange meson channel 2y , the number and the symmetrization in terms of the 2y interpolators are given for the light meson 



channels 2j, and 2j, 
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Interpolator 
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TABLE VII: Pseudotensor interpolators from irrep E. The unsymmetrized interpolator, the corresponding number in the 
strange meson channel 2^, the number and the symmetrization in terms of the 2^ interpolators are given for the light meson 
channels 2^~ and 2^"*". 
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TABLE VIII: Tensor interpolators from irrep T2. The unsymmetrized interpolator, the corresponding number in the strange 
meson channel 2^ , the number and the symmetrization in terms of the 2^ interpolators are given for the light meson channel 

2Ta • 
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TABLE IX: Tensor interpolators from irrep E. The un- 
symmetrized interpolator, the corresponding number in the 
strange meson channel 2^, the number and the symmetriza- 
tion in terms of the 2^ interpolators are given for the light 



meson channels 2^+ 



2. Baryon interpolators 

The baryon interpolators are slightly more compli- 
cated, the construction in each channel is given by: 
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Subsequent numbering of the interpolators with respect 
to gamma and smearing structure is performed in each 
channel, leading to the interpolator numbers given in Ta- 
bles |X] and |Xll Projection to definite parity is performed 
by the projection operator P± = (1 ± 7t)/2. C is the 
charge conjugation operator, in the chiral representation 
it can be written as C = 17274. The A and 17 interpola- 
tors are projected to spin | an d averaged over the three 
spatial vector components j44| . 
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TABLE X: Baryon interpolators for nucleon, E and H chan- 
nels. The Dirac structures, the quark smearings and the cor- 
responding interpolator numbers are given, x labels the 
Dirac structure of the baryon interpolators. In the nucleon 
channel we prune interpolators which are very similar to oth- 
ers, obtaining six per Dirac structure, thus a total of 18 in- 
terpolators. 
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TABLE XI: Baryon interpolators for A and il channels. The 
quark smearings and the corresponding interpolator numbers 
are given. We prune interpolators which are very similar to 
others, obtaining six interpolators. 
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